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Abstract 

The predictions of the large mass hmit for the radiative decays of the known p-wave 
charmed mesons are analyzed. Special attention is devoted to the problem of gauge 
invariance of the transition matrix elements. The width ratios arising from heavy 
quark symmetry are given for the different multipole components. Finally, some es- 
timates for the rates of these decays are given, using the constituent quark model 
of Isgur, Scora, Grinstein and Wise. 
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1 .Introduction 



The strong, weak and electromagnetic interactions of mesons containing one heavy 
quark can be described in terms of an effective theory which is simultaneously in- 
variant under heavy quark spin-flavour symmetry, the chiral SUl{2) x SUii{2) and 
the [/e.m.(l) gauge group. This effective theory has recently been applied [|l], ^ to 
predict the B* B'j decay widths from the CLEO data0 on branching ratios for 
D*~^ — i> D'-){'k). There is, as yet, no similar data on radiative decays of the known 
p-wave D mesons, although the mode 1^*^(2536) -0*^7 is reported |Q as having 
possibly been seen. However, the situation might change soon if enough statistics 
becomes available. It is therefore of some interest to see what information can be 
obtained about these decays by using solely the symmetries of the large mass limit. 
In the present paper we study the radiative decays of the members of the = I"*" 
multiplet of charmed mesons to the respective charmed s-wave states. Previous work 
on the same problem has been done in [0, [^l- These authors employed a specific 
dual model and tensor meson dominance and/or made use of the SU{4:) flavour 
symmetry to extract the necessary couplings. Such an approach is clearly not com- 
patible with the large mass effective theories referred to above. We prefer instead 
to describe the electromagnetic interactions of these mesons using a method which 
explicitly displays the new spin and flavour symmetries appearing in the heavy mass 
limit. We show that the heavy quark contribution to the matrix elements of these 
transitions is, up to 0{l/mc), completely determined in terms of the same Isgur- 
Wise function ^3/2 (t^ ■ v') which describes the semileptonic decays of the B mesons 

to excited sj' = |^ p-wave D-mesonsP]Q. Ratios of partial widths are given for 
the different multipole components. Finally, the constituent quark model of Isgur, 
Scora, Grinstein and Wise ^ is used in order to make some estimates for the widths 
of these decays. The result is that the decays of the mesons are dominated by the 
magnetic quadrupole M2 mode while the mesons decay predominantly through 
an electric dipole El mode. 



2. The method 

We are interested in the matrix elements of the electromagnetic current 

J^.m. = l^l^U - id-f^'d - ^SJ^S + '^-n^C (1) 



taken between one of the s-wave charm meson states and one of the two degenerate 

3- 
2 

M% = {M,[v')\J^,,^\Mr[v)), (2) 



states of the sY = § p-wave multiplet: 



^More precisely, this is only true for the El component of the transition. 
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where Mi{v') , M** (v) are generic notations for the respective meson states, 
(Mi,M2) = {D,D*) and (M**,M**) = {Di,D*). Their fields can be combined, as 
usual |T0[ into a Dirac matrix^ 

H = i^[D*V-/^75] (3) 
= ^{/^r'7A-^/^h5[3^7.--7A(7"-^")]}. (4) 

We will choose to evaluate the matrix element (2) of the first three terms in (1) 
(the light quarks' electromagnetic current) in a somewhat different way from the 
last one (the heavy quark electromagnetic current). Namely, the former is given 
in the effective theory by the matrix elements of the most general operator which 
includes the fields (3-4), has positive parity and is gauge invariant. To second order 
in the photon momentum it can be written as 



n 



,Jt^,Ht = fiTT[HT-^PF^p] + f,TT[Hn^F^p] (5) 
+ j^Ti[HT''dxa^^F^p] + |^Tr[i5T"9„a^,F^'^] + ^Tr[iyT"(7^t;^ + j^v^)dxF'"'] . 

XXX 

Here F'^'^ = d'^A'^—d'^A'^ is the electromagnetic field tensor and is the polarization 
vector of the emitted photon. Note that the velocities of the two heavy meson fields 
in (5) must be taken to be equal v = v', since otherwise the invariance under the 
SUspin{'2) heavy quark spin symmetry group is lost. Terms with one derivative acting 
on the heavy meson fields are forbidden by reparametrization invariance constraints, 
as shown in for the case of the strong interactions of the same mesons. /, h 
and j are unknown constants which have to be determined from a comparison with 
experiment. A^=l GeV suppresses terms with one derivative on the photon field. 
The matrix element (2) of the last term in (1) will be evaluated by using usual 



HQET methods [T^. The vector current cy^c is matched in the HQET to 0{l/m^ 
onto the operator 

JLa.y = e'-^^^'-^^ih^^^h^:^ - l^h}$[P Y - Y mi""^} . (6) 

The two velocities v, v' are taken to be different here and satisfy 

^P3/2^M = Msv'^ + K (7) 
with the photon momentum. Defining the binding energies 

K' = Ms-m,, A = Mp3/^ - m, , (8) 



factor of ^/Ms{^Mp^^^) is absorbed into the field H (T"). 
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we have 



/T Tnr A-A' A(A-A') ,^,M^ A'\, 
^ ^, _ (A-A')2 ^^^A3 

It can be seen from (10) that v ■ v' differs from 1 only through terms of 0{l/m1). 
This provides the justification for our assumption about the equahty of the two 
velocities v and v' built into the equation (5). However, in considering the matrix 
elements of the heavy quark current (6), we will stick to a two-velocities description, 
as a means for introducing the photon momentum via (7). 
Using the usual trace formalism, we have 

mv')\Jl,^y{x)\M;*{v))= (11) 

where only the leading term in (6) has been retained. Mi and M**° are the usual 
interpolating fields associated with the respective meson fields: 

{Ml , Ml) = (75 , ^*), (12) 

(Mr , Mr) = (^--^e^5[3^7," - 7a(7" - v^] , • (13) 

Unfortunately, the matrix element (11) is not gauge invariant. Contracting equ. (11) 
with = Mp^^^v^ — Msv'^ gives a result of order 0(m°), proportional to (A — A'). 
It is not difficult to trace the origin of this difficulty: it is the nonconservation of 
the HQET vector current (6) when only the first term is retained. In order to have 
vector current conservation, the second term in (6) must be included as well. Only 
then we have 

d,JLa.y = Oil/m,) . (14) 

Of course, vector current conservation holds true up to any given order in l/rric, 
provided the necessary additional terms are added to (6). 

On the other hand, the inclusion of the second term in (6) requires also consid- 
eration of the ^ — terms in the heavy quark effective lagrangian 

Chqet = h'^Kiv ■ D)h^:^ + ±-U:\iDfh^^ - -£-U:^a ■ Chi^^ , (15) 

which could give in principle contributions of comparable magnitude. As is well 
known since Luke's paper|13[ all these bring along new, unknown form-factors. We 



3 



will show, however, that to first order in —, no new form-factors show up besides 

We shall first evaluate the contribution of the second term in equ.(6). To this 
end consider the matrix elements of the operator —ih!"^) Dp ^h^\ which are given 
by 

^MsMp^^^ Tr [M/ l±^r i±^M;- {hg^, + hv'^v^ + hv'J^ + /4<7p)] (16) 

Here /i_4 are functions of v • v' . One gets a first constraint on these unknown 
form-factors by contracting with v''' and using the equation of motion h^^) D -v' — 
0: 

fi + f2V ■v' + h-U = 0. (17) 
A second constraint is obtained by taking the derivative 

= (AV - Kv)p{M,{v')\h^)Th^^\M**{v)) 

= {M,{v')\ - ih^) Dp r/iW - ih^)T Dp h^^\M**{v)) . (18) 

Contracting with v^^ employing the equation of motion v ■ Dh^^ — and comparing 
with (16), a second constraint is obtained 

/2 + fsV -V'-U^ [P^V ■ V' - mi/2{v ■ V') . (19) 

From (17) and (19) we get 

A(1) = -(A'-A)Q,(1). (20) 

The matrix element of the second term in (6) can be expressed with the help of equ. 
(18) as 

-i{M,{v')\h!~^(p 7^ - 7, P)h^^\M*;{v)) = (21) 

-2i{M,{v')\h^) Dp h^^M**{v)) - (AV - Kv)p{M,{v')\h^)-f,Yh^^M**{v)) . 

The first term on the r.h.s. is given by a formula like (16) with F = 1, where it can 
be seen that for Vp = v'^ only the term with fi survives because of the orthogonality 

relation M**°'Va = 0. But from (20) one sees that a.t v ■ v' = 1, fi can be expressed 
only in terms of ^3/2(1)- This is essentially the explanation for our simple final result. 
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Now for the contribution of the non-leading terms in the HQET Lagrangian 
(15). The h^^^\iDYh^fi term just " renormahzes" ,^3/2 by effectively changing 

il/2{v ■ v') ^ il/2{v ■ v') + —^{v ■ v') (22) 

and the chromomagnetic term h^^^a-Gh^^ gives a contribution to the matrix element 
of interest equal to 

•m,--'--' 2 2 2 

+ ^MsMp^^^ i-Tr [Mt i^^pA 1±1L^^ ^M;*°< {xicr.x - iX2{vpix - ^;a7p))] -(23) 

7], Xi and X2 are real unknown functions. There are thus five unknown form factors 
besides the Isgur-Wise function ^3^2 which determine the matrix elements of the 
(approximatively) conserved vector current C7^c in the heavy quark effective theory 
up to (9(l/mc)0. In the expressions for these matrix elements we insert v'^ as a 
function of and from (7), the meson masses Ms and Mp^^^ are expressed in 
terms of A and A' and the result is expanded in powers of In this expansion, 
one power of k^ counts as much as one power of A , A', /i_4, r/ or xi,2 as they are of 
the same order of magnitude. The leading terms of this expansion look as follows: 

{D*{v\e,)\Jl^'^^y\Dl{v,e,)) 

2 . . A -A' 



iMsMp^^^ ( -^efel^k^v, + 2— — ereL ^3/2(1) , (24) 
{D*iv',e2)\Ji^'^^'mv,e^)) = ^ MsMp^^^^ie^p,,e^,efk'<C,/2{^) , (25) 
{D{v')\Jl^^^y\Dl{v, ei)) = ^MsMp^^^ ^^e„,,,ef ^"^"^3/2(1) , (26) 



It is apparent that all the matrix elements (24-27) can be expressed only in terms of 
^3/2(1) and that none of the five subleading form factors contributes. An important 
role is played here by equ. (10) which prevents the appearance of terms proportional 
to the derivative of ^3/2(2^) at x = 1. These matrix elements exhibit approximate 
gauge invariance (up to 0{l/mc)), as shown in (14). In practice this is slightly 
inconvenient, so that we will henceforth replace A— A' by f -A; in the above equations, 
which only changes the result by a next-to- leading quantity, as is apparent from (9). 
After this change our matrix elements are explicitly gauge invariant. 

•^Similar results were obtained also in the Ref. ilql. 
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3. Results 



The four possible electromagnetic decays of the sY = | charmed mesons to the 
= I charmed mesons have, in the infinite mass limit, the following multipole 

content: D2 D*, ^ D, Di D*, Di D. 

The D2 D* transition can in general proceed also through a E3 mode, as 
far as angular momentum and parity are concerned. In the infinite mass limit 
this is forbidden. However, this prediction is not very likely to be easily tested 
experimentally. 

More interesting are the heavy quark symmetry predictions for the ratios of 
partial amplitudes for these decays: 

Aei{D; ^ D*7) : Aei{D; ^ D7) : Aei{Di ^ D*^) : Aei{D, ^ D^) 

= V3 : : 1 : V2 (28) 

AM2iD; -> D*7) : AM2iD; ^ D7) : ^^2(1^1 ^ D*-f) : Am2(/^i -> D-f) 

= V3 : V2 : V5 : 0, (29) 



which can e.g. be obtained along the lines of Ref.|]T4[. Using (5) and the matrix 
elements (24-27) we obtain the following one-photon widths, which automatically 
satisfy the above relations: 

T[Di^D-i) = -^ji—\k\{eqF + eQ^ ) + a——\k\ [egG) , (32) 

T{D,^D^) = + ■ (33) 

Here and eg (= | for Q = c) are the light and the heavy quark electric charges 
and F , G are given by 

F = fi-f2 + ^{gi + 2g2 + 4g,), (34) 

X 

G = -^{g, + 2g2). (35) 

The heavy quark contribution to the M2 partial width has been written only for 
the transition D2 —>■ -D7, since it is only in this case that it can be expressed solely 
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in terms of ^3/2(1)- Fo^^ ^ill the other decays, it depends also on some of the five 
unknown subleading form factors. 

We will use in the following the ISGW valued] 

e3>(l) = 0.584. (36) 

This corresponds to a transition c ^ c and is slightly different from the correspond- 
ing value for a 6 — c transition, of 0.537 A similar value is also obtained from 
a recent QCD sum rule calculation | p. 5||. Although this value is the result of a model 
calculation, 1^3/2(1) can be obtained experimentally, from a study of the differential 
cross-section for the semileptonic decays of B mesons into p-wave D mesons. In this 
sense the relations to be derived below are truly model-independent. 

One first set of predictions is nothing else but the mass corrections to the ratios 
(28-29). In the approximation that F and G are constants independent of the photon 
energy, these are: 

TEiiDlt -> D*/^) : TEiiDl^ - Dtl) : r^i(/^+ - Z^^t) : ^EiiDt, ^ D+7) 

= 3.694 : : 1 : 4.071 (37) 

rM2{D:t - D:^^) : TM2{D:t - D^^) : TM2{Dj, ^ Dl^^) : 1^2 (/^+i - Df^) 

= 0.965 : 2 : 1.212 : 0. (38) 

and similarly for the mesons: 

= 3.539 : : 1 : 4.022 (39) 
Tm2{D;'' ^ D*'^) : Tm2{D;'' - D^^) : Tm2{D^i ^ D*'^) : TM2{D'i ^ D^) 
= 0.964 : 2 : 1.125 : 0. (40) 

In the corresponding one pion decay case a small mixing between the Di and the 
axial vector member of the sj' = |^ multiplet could significantly alter the predicted 
amplitude ratio. However, in our case this is no longer true. Both these states can 
decay to the s-wave mesons through a El mode, and as far as the M2 channel is 
concerned, it is only the decay of the vector member of the sj^ = multiplet which 
is forbidden by large mass limit selection rules. 

It is clear that without any knowledge of the constants F and G no further ad- 
vance can be made. Since the value (36) has been obtained by using the constituent 



*The relation to the function f^^2 used in g is = V3'?3/2 



^Notc that what we call S,^^ the usual Isgur-Wise function, but rather its renormalization 

group invariant version, which includes a scaling factor with the respective velocity dependent 
anomalous dimension and has a logarithmic dependence on the heavy quark mass. This in part 
explains the two different numerical values quoted above. 
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quark model of Isgur, Scora, Grinstein and Wise, consistency requires that we use 
the same model for determining them. Although the authors of Ref.0 made use 
of their model only to calculate matrix elements of heavy quark current bilinears, 
it can give matrix elements of light quark currents as well, thereby providing a way 
for evaluating our parameters F and G. The model is well tested and is known to 
give reliable results for both the s- and p-wave charmed mesons. Furthermore, it 
provides a simple way for including SU{'i) violating effects induced by a constituent 
strange quark mass nts different from m„, m^. 

To determine F and G for the meson case, consider the following matrix 
element of the light quark current uj^u, which according to the effective lagrangian 
(5) is given by 

{D'{v)\u^^u\D',iv,e)) = A^MDoMj,o[-iv ■ k)e^ + (e ■ A;)t;^]FBO . (41) 

The same matrix element is written as 

{D^{v)\u-i^u\D^^{v, e)) = ^MdoMdo[-2.673{v ■ k)e^ + 2.396(e ■ k)v^] (GeV) , (42) 

where we make use of the relations given in the Appendix B of Ref.0. Here some 
care is needed, as the light quark in a D meson is an antiquark. Therefore the direct 
application of the relations of Ref. gives rather the matrix element between the 
corresponding D mesons. However, these can be related by crossing to the matrix 
elements of interest. Unfortunately, the transition matrix element (42) is not gauge 
invariant. This is a typical problem of constituent quark model calculations of 
radiative transition matrix elements (see e.g. |]I^) and we will deal with it by 
adopting as the value of the coupling 4F the average of the two numbers in (42), 
with an error given by their difference. This should give at the same time some 
measure of the model dependence of the result. Thus we take 

Fdo = (0.634 ± 0.034) Gel/"^ (43) 

In obtaining this number we have used a value k = 1 for the "relativistic correction 
factor" K (see Ref.23 in 0). Adopting k = 0.7 lowers F by a factor of about 0.5. 
Similarly the value of G can be obtained by calculating the matrix element 

{D\v)\u^^u\D;\v, e)) = -2t^MnoMj,,oGDoe^,xpe''''k^v^k' (44) 

with the result 

Gdo = -2.1Q8 GeV-\ (45) 

Inserting these values for F]jo and G^o in the rate formulae (30-33), we obtain the 
partial widths in Table 1. Here a charmed quark mass rric = 1.82 GeV has been 
used. 
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A similar calculation gives the following values for the corresponding couplings 



for the mesons: 

= (0.493 ± 0.083) GeV'^ (46) 
= -1.254 Ge^-^ (47) 

In obtaining these numbers the following values for the oscillator strength parameters 
of the model have been used: Ps = 0.40 GeV and j3p = 0.35 GeV. For the flavour 
content sc these parameters are not available in the Ref . and therefore we simply 
took them slightly larger than for uc, as seems to be the case for the light mesons. 



The mass of the tensor state has been taken equal to M^*+ = 2564 MeV |T8 
The partial rates for this case are also shown in the Table 1. Several remarks are in 
order about these results: 

• The heavy quark contribution amounts to about 34% in absolute value in the 
El amplitude for the mesons' case and 56% for the case. In the former 
it has an enhancement effect while in the latter case it contributes with an 
opposite sign. On the other hand, in the M2 amplitude for the D2 D 
transition the heavy quark current contributes negligibly (under 4%). Very 
likely, the same is true for all the other M2 amplitudes. 

• Our estimates for the total one-photon widths agree well with the predictions 
in for the decay D*~^ and are lower for the same decay of the 
meson. Compared to the our estimates for T{D2 — > -D7) are down by a 
factor of 1.8 and 5.5 for the and cases respectively. Also, we get a 
value for r{D2^ —>■ D*^'~f) which is lower by a factor of 5.2 than the one given 
in 1^. As a general common point we mention the strong suppression of the 
decay widths of the mesons compared to the mesons. 

• These estimates yield branching ratios of about (0.9-1.5)% for the radiative 
decays of the mesons and larger than (0.02-0.4)% for the Dfi, for which 



an upper limit for the total width of 3.9 MeV exist s [[l^] . 



The largest source of errors in our approach seems to be the neglect of higher 
order terms in the effective lagrangian (5), with more derivatives on the photon 
field. Such terms would be a source of additional photon-energy dependence of the 
couplings F and G. An estimate for this dependence can be obtained by formally 
setting = in (34). Then the dependence of F with the photon energy can be 
expressed simply in terms of G. This yields a variation of F among the three allowed 
El transitions of about 15 %, with the attendant corrections to the amplitude ratios 
(37) and (39). 

The finite mass corrections are only significant in the El amplitudes, and we 
have estimated them to contribute up to 10% of the amplitude. Much more likely 
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are important the corrections due to SU{3) violation. These can be calculated in 
a systematical way|l|] using the heavy hadron chiral perturbation theory. We have 
only evaluated them in an effective way, from a constituent quark model, to be of 
about 25% for F and 50% for G. 

Similar methods can be applied to the description of the electromagnetic decays 
of the members of the sj' = |^ multiplet of charmed mesons and of the similar 
excited B mesons. Since the former are expected to be quite broad (F ~ 200 MeV) 
because of the pionic decay mode, their photonic branching ratia are surely very 
small. As for the latter, in the view of the absence of any experimental data about 
these states, such applications must be reserved for the future. 
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Decay 


r^;i (keV) 


(keV) 


^tot (keV) 


1^1 (MeV) 


r(DO ^ 

T{Df D^-f) 


(245 ± 18) 

(60 ± 5) 
(222 ± 16) 





101 

87 
181 


(245 ± 18) 
(161 ±5) 
(309 ± 16) 
181 


494.93 
381.05 
410.38 
522.63 


r(D+ ^ d:+7) 
r(L>:+ ^ L>r7) 
r(D:2+ - i^,+7) 


(1.6 ±2.3) 
(0.4 ± 1.0) 
(1.4 ±2.0) 




10.0 
8.0 
16.0 


(1.6 ±2.3) 
(10.4 ±1.0) 
(9.4 ±2.0) 
16.0 


503.77 
389.97 
413.56 
526.12 



Table l.The calculated partial widths for the radiative decays of the sj^ — | 
p-wave charmed mesons to s-wave charmed mesons. In the last column the corre- 
sponding photon momentum is given. 
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